Electron-electron interaction effects in quantum point contacts 
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We consider interaction effects in quantum point contacts on the first quantization plateau, tak- 
ing into account all non momentum-conserving processes. We compute low-temperature linear and 
non-linear conductance, shot noise, and thermopower by perturbation theory, and show that they 
are consistent with experimental observations on the so-called "0.7 anomaly". The full temperature- 
dependent conductance is obtained from self-consistent second-order perturbation theory and ap- 
proaches ^ e^ /h at higher temperatures, but still smaller than the Fermi temperature. 

PACS numbers: 72.10.-d, 73.23.-b, 72.10.Fk 



CD 



o 



> 

ON 
OO 
ON 

o 
1^ 
o 



X 



Conductance quantization in a quantum point contact 
(QPC), first observed in 1988 [l|, constitutes a classic 
textbook effect of mesoscopic physics. On top of the inte- 
ger conductance plateaus G = nGo (where Gq = 2e^ /h) 
observed as a function of gate voltage Vg, many exper- 
iments have pointed to the existence of a so-called "0.7 
anomaly" in the conductance and other transport quan- 
tities [l, H, 0, H, 111 ■ Most prominently, the 0.7 anomaly 
implies a shoulder-like feature in the conductance G{Vg) 
around G k, 0.7 Gq seen at elevated temperature T (or fi- 
nite voltage V) near the first quantized plateau 0, H, 01 j 
accompanied by a shot noise reduction Q. Given the 
conceptual simplicity of a QPC and the fact that the 
0.7 anomaly has been observed in a variety of material 
systems by different groups over more than a decade, 
it is quite amazing that still no generally accepted mi- 
croscopic theory exists, apart from an overall consensus 
that one is dealing with some spin-related many-body 
effect. Such a theory should be able to explain all the 
experimental data in a unified and physically consistent 
manner. 

While phenomenological models Q , assuming the exis- 
tence of a density-dependent spin gap, can provide rather 
good fits to experimental data, the presumed static spin 
polarization due to interactions within the iocaJ QPC re- 
gion is not expected in the presence of unpolarized hulk 
reservoirs. Recently it was also pointed out that spin 
symmetry-broken mean-field theory is unable to recover 
the correct T dependence of the conductance 0]. Other 
proposals assume the existence of a quasi-bound state in 
the QPC region, leading to a Kondo-type scenario as en- 
countered in transport through interacting quantum dots 
[§, Such a quasi-bound state was indeed found in 
spin density functional theory (SDFT) calculations Q, 
but other SDFT works did not reach such conclusions 
pd| . Further proposals involve phonon effects [l^- Sev- 
eral publications have suggested that taking into account 
only electron-electron (e-e) interactions may result in a 
reduced conductance at elevated temperatures, without 
the need for additional assumptions of spin polarization 
or a localized state [l^, 14, 15, 16, JJi]. However, a physi- 
cally consistent picture explaining the temperature, volt- 
age, and magnetic field dependence of the conductance. 





FIG. 1: (Color online) Illustration of the two-electron non 
momentum-conserving scattering processes that give rise to a 
correction to the transport properties at the beginning of the 
first plateau. The full (black) lines represent incoming elec- 
trons, while the dashed (red) lines are the outgoing electrons. 
The thick (blue) lines define the edge of the QPC. Only scat- 
tering between different spins is present to leading order in 
T/Tf due to the Pauli principle. 



as well as thermopower and noise experiments, is still 
lacking. In this paper, we show that a careful consider- 
ation of non momentum-conserving e-e interaction pro- 
cesses in QPCs may allow for a consistent theory of the 
0.7 anomaly. 

The lack of momentum conservation in e-e scattering 
processes is due to an emerging lack of translational in- 
variance relevant for the low density regime kpL ^ 1, 
where kp is the local Fermi momentum and L a typical 
interaction length-scale of the QPC (see below). There- 
fore this effect is dominant at the onset of the plateau 
and gradually disappears for larger electron density in 
the QPC. We will focus on the first conductance plateau, 
where the QPC has only one open channel (ID mode) 
[isl . Iigt . Now e-e interactions give the contribution 

<t<t' 

(1) 

to the Hamiltonian, where '^a{x) is the ID electron field 
operator for spin a =t,i- The pair potential W{x,x') 
takes into account screening processes due to closed chan- 
nels and nearby gates, as well as semiclassical slowing 
down [l7j , and therefore depends not only on the rel- 
ative coordinate x — x' but also on the center of mass 
X ^ {x + x')/2. In fact, the range of the interaction 
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W{X) is determined by the above length-scale L. The 
e-e interaction now enters the ID description in two dif- 
ferent ways. There are (i) one-particle effects described 
through a self-consistent potential, given by the real part 
of the self energy (in the simplest case, this is the Hartree- 
Fock potential). At T = 0, the imaginary part of the self 
energy for our Fermi liquid starting point is zero [20| . 
and the plateau then occurs where the open channel suf- 
fers no backscattering. This part is thus already con- 
tained in the potential forming the QPC. In addition, at 
finite T, we have (ii) inelastic e-e scattering processes, 
which are the focus of our work and give rise to a re- 
duced conductance by changing the number of left- and 
right-movers, and hence the current. This is illustrated 
in Fig. [U where two important types of non momentum- 
conserving scattering events are illustrated. Process b(2) 
describes the simultaneous backscattering of two elec- 
trons with opposite spins and has been discussed on a 
perturbative level in Ref. [3]. The process b(l), where a 
single electron is backscattered, has not been studied be- 
fore. On top of these two, there are e-e forward scattering 
and backscattering processes (momentum-conserving in 
a long ID wire). While high-T transport properties are 
affected by all interactions, we show now that the leading 
low-T behavior is fully determined by the two processes 
in Fig. m 

Let us calculate the non-linear conductance, the ther- 
mopower, and the shot noise to leading order in the in- 
teraction. Perturbation theory gives the interaction cor- 
rection to the current (here, h = 1 and e > 0) 
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= 1- (Afc(i)+A(2))(7rT/T;^) 



(2) 



(A(i)/4 + ^(2)) {eV/epf + 0(W^'), 



M1.2) = W^6(i,2)^F/(487r24) 



with coefficients A, 
responding to the b(l) and b(2) processes, where 
ep = ksTp is the Fermi energy and M^(,(i.2) = 
J dxdx' iy(a;^a;')e'fcf'(=^+^')eifei^(^T^') . Already at this 
point, we observe a correspondence to experimental ob- 
servations, namely the reduced conductance I /V with 
increasing T and/or V . Furthermore, to leading order in 
T/Tp, where all scattering happens at the Fermi level, 
only opposite spins interact due to the Pauli principle: 
for equal spins, the exchange term tends to cancel the 
direct term. As a consequence, we can also understand 
the behavior at large magnetic fields, where the T = 
plateau occurs at e^ /h. In that case, to order (T/Tp)'^ 
no interaction renormalization of the conductance arises 
PH . This is consistent with experiments, where no sup- 
pression is observed at the half-plateaus. 

Another experimental observable probing the en- 
hanced phase space for e-e scattering at higher T is the 
thermopower S{T) p^ . for which perturbation theory 
predicts 
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Since the non-interacting thermopower is exponentially 
small [cx exp(— Tp/T)] at the conductance plateau, the 
interaction correction completely determines the low- 
temperature thermopower [2^. The enhanced ther- 
mopower (as compared to the non-interacting one) is in 
qualitative agreement with experiments at the anomalous 
plateau 0]. 

Next we calculate consequences for another observable, 
namely non-equilibrium noise. The zero-frequency shot 
noise follows from the (symmetrized) two-point correla- 
tion function of the current operator. Perturbation the- 
ory yields for the backscattering noise power 

SBiV,T) ^ 2e [24,(2) (l/,r)coth(el//fcBr) (4) 

+ /6.(i)(y,T)coth(e1//2fcBT)], 

where Ibs(i.2) £^re the current corrections due to Wft(i 2) 
quoted in Eq. ^ (defined positive for V > Q). This 
is nothing but the famous Schottky shot noise relation, 
encoding the charge of the backscattered particles. Equa- 
tion ^ predicts an additional factor of two for the 
b(2) contribution, because two electrons are backscat- 
tered in that event Direct calculation then yields 
the full noise power of the transmitted current as St ~ 
Sb + AGokBT - SkBTdvhs; where = 4,(1) -I- 4,(2) . 
Recent noise measurements on the first quantized plateau 
were compared to the corresponding single-particle pic- 
ture 0, and a reduced noise power was observed on the 
conductance anomaly. For that comparison, one sub- 
tracts the thermal noise and defines the excess noise as 
Si ^ St ~ 4:G{V,T)kBT. For a non-interacting system, 
SfP = 2GQR{eV coth{eV/2kBT) - 2A:bT} to lowest or- 
der in the reflection coefficient R = Ibs/GoV, see Ref. 0. 
Thus the difference between the true excess noise and its 
single-particle value is 



Si - SfP 
2GoeV{T/TFy 



-2A 



6(1) 



eV 
k^' 



-A(2) /i(el//fcsT), (5) 



where /i(.t) = — 8a;+ (tt^ +a;^) tanh(a;/2). This expression 
shows that for eV < 6.507 fcsT, regardless of Af,(i 2), 
the measured noise is always smaller than predicted by 
a single-particle analysis. This situation corresponds to 
the experimental work of Ref. @, where eV ^ Sfc^T. 

It is clear from all these perturbative results that for 
low energies, 1/, T — > 0, all interaction effects disappear. 
The perturbation theory results presented above, how- 
ever, obviously break down at higher temperatures or 
voltages. From Eq. we find the temperature scale for 
this crossover to a strong-interaction regime. 



kBT* 



cx 



Wbkl ' 



where A 
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v4f,(2) and 



(6) 



Con- 



(A(i)+A(2))(T/Tp)' 



(3) 



trary to the usual situation encountered in mesoscopic 
physics, the nontrivial question to be answered thus con- 
cerns the high-temperature limit (but still T <^ Tp). To 
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make progress in the relevant temperature regime 



(7) 



let us consider a simplified model pair potential entering 
Eq. 11]), see Ref. [13, 



Wix,x') =VoSix)Six'), 



(8) 



which implies Wi,{i) = Wb(2)- We express the interaction 
strength in the dimensionless parameter A = 
Estimates for A in GaAs heterostructures gives Awl, sec 
Rcfs. Ii,[l3, which then yields keT* /Tp « 0.1 allowing 
for a study of the temperature range ([7]). 

In order to treat the local interaction ([8]), we start from 
the Dyson equation for the full Keldysh single-particle 
Green's function (OF) G{x,x' ]uj), 

G{x, x'] Lo) = Ga{x, x'\ w) + Go(x, 0; w)S(tj)G(0, x'; w), 

(9) 

which is a 2 X 2 matrix in Keldysh space. The self energy 
due to Eq. ([5]) acts only at a; = x' = 0. Note that = 
Gj^, i.e.; wc can suppress the spin index. For simplic- 
ity, we now assume a parabolic dispersion, = fc^/2m, 
for the open channel. The charge current operator is 

^ = ^E. [^U^)d.'^A^) - id.H(.^))^A^)], and we 
evaluate (/) at x = 0, where it can be expressed in terms 
of the local GF G(a;) = G(0,0;w) and the self energy 
S(ti;). In fact, some algebra shows that only the local 
spectral function A{u!) = —2 lmG^{ui) enters the current 
formula for the contact interaction (ISl, 



/ = 



2e 



^oiw) 



(10) 



where /^/^ are Fermi functions in the right/left lead, 
and Aq{ll}) = 2'Kd{uj) is the non-interacting spectral func- 
tion. Here, d{Lo) is the density of states 2'Kd{uj) ~ 
{2ra/uj)^/'^9{u}). Remarkably, the nonequilibrium current 
through the interacting QPC is thereby fully expressed 
in terms of the local retarded GF only. So far, the given 
relations are exact, but to make progress, one needs to 
approximate the self energy. We take the full second- 
order self energy^ 

E'-(a;) = dt e"^'[G<{~t)G>{t)G>{t) 

-G>i-t)G<it)G<it))], (11) 

and make it self-consistent by using the interacting 
(Icsscr/greatcr) GFs. The corresponding diagrams are 
shown in the inset to Fig. [21 The approximation (|lip is 
the simplest way to describe equilibration between left- 
and right-moving electrons in an interacting QPC. In 
what follows, we confine ourselves to the linear conduc- 
tance regime, where the spectral function in Eq. (jlO[) 
can be calculated in equilibrium by solving Eq. pT|) . 
and where we can replace /j? — /" eV[—dujf{uj)], 
where f{uj) is the Fermi function. In linear response, the 



lesser/greater GFs can be written in terms of the local 
spectral function A(uj), 



G</>{t) = ±1 



duj 
2^ 



*A(c.)/(±t.) 



(12) 



This suggests a natural iterative way to self-consistently 
solve for the conductance: Starting with the initial guess 
A{u!) = Aq{lli), one computes S''(w) from Eq. (fTTI) . which 
in turn defines a new retarded GF and a new guess for 
A{ijj). This procedure is iterated until convergence has 
been reached. For the parameters below, this numeri- 
cal scheme is convergent and can be implemented in an 
efficient manner. 

The numerical results, for A = 0.8 shown in Fig. [21 
accurately reproduce the above perturbative results at 
low r, but also allow to cover the interesting high- 
temperature limit. Our data for different A fall to high 
accuracy on the simple function 



G{T) 

Gq 



= b 



1 - b 



(13) 



where b sets the high-temperature saturation value. 
While this functional dependence is somewhat similar 
to the phenomenological Kondo-type function used in 
Ref. [3, our numerical data fit better to Eq. (fT5)) . It is 
also possible to obtain equally good fits to the activated 
T dependence reported in Ref. Q , see also Ref. 0] , 



G{T) 
Go 



1 - (1 - a)e 



(14) 



where a again denotes the high-T limit. The values for 
and extracted from best fits to our numerical data 
are summarized in the inset in Fig. [21 Remarkably, both 
temperature scales are of the same order. Moreover, they 
are lowered by increasing the interaction strength A. For 
high r, the conductance appears to approach the satu- 
ration value G ~ /h. Similar saturation value has also 
been reported for long wires [l3j . with the same T = 
conductance Go- The new feature for QPCs comes from 
the non-momentum conserving interactions, resulting in 
a distinct low-to-intermediate temperature dependence 
G{T). The perturbative correction is not present in 
the long wire results [ll| , but is seen experimentally [J . 

As a final remark on the numerical solution of the self- 
consistent approach, wc mention that thcrmopower (data 
not shown) exhibits a crossover from the S oc law at 
low T, see Eq. ([3]), to a linear- in-T behavior at elevated 
temperatures. 

It is also instructive to discuss our model in terms of 
an Anderson model. For the model pair potential ([H]), 
by spatial discretization our Hamiltonian maps to a ID 
tight-binding chain with hopping matrix elements t and 
on-site interaction U acting at one site {x = 0) only [2^ . 
We thus arrive at an Anderson-type impurity model, sim- 
ilar to the one used in Ref. [9| to describe interactions in a 
QPC in the Kondo regime. However, we consider a rather 
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0.1 0.2 0.3 0.4 

T/Tf 

FIG. 2: (Color online) Temperature dependence of the lin- 
ear conductance G for A = 0.8. Dots denote self-consistent 
numerical results, the solid curve gives the perturbative esti- 
mate ([21, and the dotted curve is a fit to Eq. (|13p . Inset: Fit 
parameters entering best fits of Eqs. (|13fl and (|14|l to numer- 
ical data for T/Tf < 0.4. Note that a and h are somewhat 
different. The lower inset shows the two second order self- 
consistent energy diagrams in Eq. (jlip . 

different parameter regime, where U is of the same order 
as the hybridization T and can be parametrically larger 
than the bandwidth D ~ \t\. Employing Eq. ([5]), with 
ep « D, the interesting temperature range ([T]) translates 
to D"^ /U <C ksT ^ D, where our claim is that G ap- 
proaches w /h. While the Kondo model requires the 
formation of a local moment, this is not the case for the 



present approach. Instead, our high-temperature limit 
may be described as an incoherent Fermi liquid, with full 
relaxation between left- and right-movers. In fact, one 
can establish that a simple Boltzmann-type approach has 
a high-tcmpcrature solution where the out-going distri- 
bution function is a mixture of the incoming left- and 
right-mover's distributions. Such an Ansatz leads to the 
conductance formula (fT5|) with b = 1/2, which is the rea- 
son for using that as a fitting formula. Unfortunately, 
a Boltzmann approach is conceptually difficult to justify 
due to an inherent normalization problem [T7|, i.e. one 
cannot define a proper local distribution function in k- 
space for this model. 

In conclusion, we have considered interaction effects 
in short QPCs and shown that taking into account non 
momentum-conserving processes, we can qualitatively 
account for the experimentally observed behavior of the 
linear and non-linear conductance, thermopower (includ- 
ing their magnetic field dependencies) and shot noise 
at the so-called 0.7 anomaly. The gate voltage depen- 
dence can also be explained within the present scheme, 
because the backscattering is suppressed for larger values 
of kpL. In the high-temperature (but still T ^ Tp) non- 
perturbative regime, our second-order self-consistent ap- 
proach predicts that the conductance approaches w e^//i. 
It is an open and interesting problem to verify this result 
by other non-perturbative methods. 
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